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Background: A component of azimuth correlations from high-energy heavy ion collisions varying 
as cos(20) and denoted by symbol vi is conventionally interpreted to represent "elliptic flow," a hy- 
drodynamic manifestation of the initial-state A-A overlap geometry. Several numerical methods are 
used to estimate V2 , resulting in various combinations of "flow" and "nonflow" that reveal systematic 
biases in the vi estimates. QCD jets contribute strongly to azimuth correlations and specifically to 
the cos(20) component. Purpose: We question the extent of jet-related ("nonflow") bias in and 
hydrodynamic "flow" interpretations of V2 measurements. Method: We introduce two-dimensional 
(2D) model fits to angular correlation data that distinguish accurately between jet-related correla- 
tion components and a nonjet azimuth quadrupole that might represent "elliptic flow" if that were 
relevant. We compare measured jet-related and "flow" '-related data systematics and determine the 
jet-related contribution to V2 measurements. Results: Jet structure does introduce substantial bias 
to conventional U2 measurements, making interpretation difficult. The nonjet quadrupole exhibits 
very simple systematics on centrality and collision energy — the two variables factorize. Within a 
Au-Au centrality interval where jets show no indication of rescattering or medium effects the nonjet 
quadrupole amplitude rises to 60% of its maximum value. Conclusions: Disagreements between 
nonjet quadrupole systematics and hydro theory expectations, the large quadrupole amplitudes ob- 
served in more-peripheral Au-Au collisions and a significant nonzero value in N-N w p-p collisions 
strongly suggest that the nonjet quadrupole does not arise from a hydrodynamic "flow" mechanism. 

PACS numbers: 13.66.Bc, 13.87.-a, 13.87.Fh, 12.38.Qk, 25.40.Ep, 25.75.-q, 25.75.Gz 
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I. INTRODUCTION 



Measurements of a Fourier component of the "az- 
imuthal anisotropy" of particle momenta in RHIC heavy 
ion collisions have been interpreted to indicate produc- 
tion of a thermalized QCD medium with low viscos- 
ity, frequently invoked as evidence for a "perfect liquid" 
( e -g-> EHUD- That conclusion is based on a conventional 
interpretation of the V2 = (cos(2 (f>)) anisotropy compo- 
nent as a measure of elliptic flow, a conjectured hydro- 
dynamic (hydro) response to density and pressure gradi- 
ents in the initial collision system corresponding to the 
transverse eccentricity of the A-A overlap region Q. In 
a hydro context large elliptic flow values combined with 
other measurements are interpreted to imply large energy 
densities, rapid thermalization and small viscosities (HHf. 

However, questions persist concerning measure- 
ments, their accuracy and their interpretation. Con- 
ventional i>2 measurements @ do not distinguish accu- 55 
rately between an isolated azimuth quadrupole (m = 2 
cylindrical multipole) Fourier component conjectured to 
represent elliptic flow and "nonflow" — a catch-all term 
representing several possible contributions to vi, but 
mainly the m — 2 azimuth Fourier component of a two- eo 
dimensional (2D) peak attributed to jets What- 
ever the precision of V2 measurements the physical phe- 
nomena actually represented by any cos(2</>) asymmetry 
measurement can be questioned [l2| LL3J ■ 



Conventional quadrupole measures W2 {method} moti- 
vated from a hydro context [3[ are difficult to interpret, 
and the statistical properties of some V2 methods lead 
to substantial systematic bias. In Ref. [l4[ it was shown 
that event-plane U2{EP} [f| is a close approximation to 
two-particle cumulant i>2{2} [l5j . in turn equivalent to 
the m = 2 Fourier coefficient of a projection of all angu- 
lar correlations onto ID azimuth [see Eq. (JTJ], which may 
include a large contribution from a prominent 2D peak 
interpreted in a perturbative QCD (pQCD) context as 
representing minimum-bias jet structure |7l. [la - flli ] . Four- 
particle cumulant t>2{4} [Uj, Hl[ reduces, but does not 
necessarily eliminate, the jet- related contribution [55j . 
Other mo del- dependent strategies have been invoked in 
attempts to reduce the "nonflow" (jet) contribution to 
v-i , but the effectiveness remains uncertain [l3|> [22[ . 

An alternative method introduced in Ref. [23[ employs 
physical-model-independent analysis to isolate geometri- 
cally the nonjet azimuth quadrupole from other contribu- 
tions. Nonjet (NJ) quadrupole amplitudes are obtained 
from fits to 2D angular correlations on azimuth <f> and 
pseudorapidity 77. Measurements of the NJ quadrupole 
over all Au-Au centralities and a large energy interval 
provide qualitatively new insights into the quadrupole 
phenomenon conventionally attributed to elliptic flow. 
NJ quadrupole amplitudes obtained with 2D model fits 
follow simple trends on centrality and energy described 
by just two initial-state parameters for all systems down 
to y/SNN = 13 GeV. 
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The quadrupole analysis method introduced in 
Ref. [23[ is based on algebraic study of v 2 methods 
in Refs. [IH [T3 and initial experience with 130 GeV 
data in Refs. [a H3| where general model-fit analysis of 
2D angular autocorrelations was first introduced. The 
same model-fit method was refined and elaborated in 55 
Ref. where the primary focus was the energy and cen- 
trality systematics of angular correlations attributed to 
minimum-bias jets or minijets. The present study com- 
bines the numerical results of Refs. @,[H1 to examine the 
systematic relation between the NJ quadrupole and mini- 60 
jets and to test the validity of the conventional "elliptic 
flow" interpretation for the former. 

In this study we examine the distinction between non- 
jet and jet-related quadrupole contributions in relation 
to other correlation structure. We review the centrality 65 
and energy dependence of the NJ quadrupole in terms of 
Glauber-model parameters as reported in Ref. [23| and 
contrast those trends with minijet systematics as estab- 
lished in Ref. [7|. We compare NJ quadrupole results 
with previous v-i {method} measurements and with hy- 70 
dro expectations. We conclude that NJ quadrupole vari- 
ations on energy above 13 GeV and all Au-Au central- 
izes are remarkably simple. Those trends and compar- 
isons with minijet systematics appear to contradict con- 
ventional hydro expectations for elliptic flow. For ex- 75 
ample, the NJ quadrupole increases to 60% of its maxi- 
mum value within a Au-Au centrality interval where the 
lowest-energy-jet-related correlations are consistent with 
a transparent collision system, as explained below. 

This article is organized as follows: Sec. [TT| reviews 
analysis methods applied to 2D angular correlations and 
"flow" analysis. Sec. [Til] presents measured angular cor- 
relations and 2D model fits. Sec. IIVI reviews model-fit 
results for jet-related and nonjet quadrupole correlation 
components and a universal parametrization of energy 80 
and centrality dependence for the latter. Sec. [V] compares 
nonjet quadrupole and jet-related trends in the context 
of hydrodynamic expectations for the former. Sec. IVII 
presents a discussion of selected results, and Sec. IVIII 
summarizes 



II. ANALYSIS METHODS 



A major emphasis of this study is accurate distinction 
between jet-related and nonjet quadrupole components 

45 of angular correlations and the energy and centrality sys- 95 
tematics of the latter — what those imply for physical in- 
terpretation of the nonjet quadrupole phenomenon. Wc 
examine the underlying assumptions and systematic un- 
certainties of the model-fit analysis method in compar- 

50 ison with alternative tj 2 analysis methods that seem toioo 
support flow interpretations. 



A. Correlation spaces 

Two-particle correlations are structures in the pair 
density on 6D momentum space (pti, f)i, <Pi,Pt2, t?2, ^2) 
that deviate from some defined reference density. In this 
analysis we study pt-integral correlations on angular sub- 
space (771, 4>\, 772, 4>2), where the angle parameters for rel- 
ativistic collisions are pseudorapidity 77 (related to polar 
angle 6) and azimuth (j). We can reduce (771, <ji>i , 772 , 4>2) to 
a viewable 2D space with no significant loss of correlation 
information by using angular autocorrelations (25|. 

An autocorrelation as conventionally defined is derived 
from a pair density p{x\,X2) by averaging along diago- 
nals in space (x\,X2) parallel to sum axis = x\ + x-i. 
The averaged pair density p(x&) on difference axis xa = 
x\ — X2 is then an autocorrelation 25j. For correlation 



structure approximately uniform on xs ( "stationarity" ) , 
typical over 2ir azimuth and within a limited pseudora- 
pidity acceptance A77, angular correlations remain undis- 
torted [24j| . Within the STAR time projection chamber 
(TPC) acceptance [26[ 2D angular autocorrelations are 
lossless projections of pt-integral two-particle momentum 
space onto subspace (?7a,<M) 0- The 4>a axis is di- 
vided into same-side (SS, \4>a\ < 7r /2) and away-side 
(AS, 7r/2 < |0 A | < tt) intervals. 



B. Correlation measures 

There are several alternatives for the definition of a 
correlation measure. The basic element is a histogram of 
covariances representing correlations of event-wise fluctu- 
ations between pairs of 2D bins on (77, <j>). Ap = p — p re f 
represents a covariance density, where object pair density 
p contains the structure of interest and reference density 
p re f may be defined in terms of a factorization assump- 
tion or constructed from mixed-event pairs. 

Per-pair density ratio Ap/p rc f = p s ib/Pmix — 1 (some- 
times referred to as a "correlation function" and denoted 
by C) varies with system size as 1/ n c h (n c h is charge mul- 
tiplicity) absent other physical changes. p s ib and p m i x 
represent sibling (same-event) pairs and mixed-event 
pairs. In previous analysis we introduced a statistical 
measure whose variation with n c h reflects only nontrivial 
physical changes in correlations, the per-particle density 
ratio Ap/y/p IC { (Pearson's normalized covariance [271 1281] 
converted to a density ratio) that exhibits the desired 
properties, since ^/p re f cx n c h 0, El E3, H^- We intro- 
duce previous correlation measurements in terms of that 
measure. We also reconsider what "particle" type best 
serves as a scaling reference in a given context — final- 
state hadrons as in Ap/ y'Pref, initial-state participant 
nucleons or the number of N-N binary collisions — and 
rescale some of the correlation data accordingly. 
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C. Two-dimensional correlation model 50 

Inspection of the 2D data histograms reveals that pt- 
integral pair-density difference A/?(?7a,0a) contains two 
types of structure: ^-dependent ID and 2D peaks and 

5 ?7A-hidependent sinusoids cos(</>a) and cos(2</>a), where 55 
the cos(20a) sinusoid (quadrupole) can be related to v 2 
measurements. We therefore define a model of 2D an- 
gular correlations that includes a part varying with 77A 
(2D) and a part independent of tja (ID) composed of the 

10 m — 1, 2 terms of a (truncated) Fourier series 60 

-% - ^(VAM+2j: ^cosWa).(I) 

VPrcf VPrcf , VPref 

v ' v ' m— 1 v ' 

No higher terms in the Fourier series are required by 
the data @, HH, [29|. Fourier coefficients Ap[m]/^/p rof = 65 
V^/27m include power- spectrum elements = 

Sr^i" 1 C0S ( TO [^i ~ <A?'D = n ( n ~ 1)(cos(to(/)a)) [3, 
15 where n is the multiplicity in one unit of 77 and 2ir az- 
imuth, so n/2ir ss d 2 n c h/dr)d(f) = pa. Angle brackets 
denote event-wise means. The first term A/^d / \fp~rej ^ s 
a combination of ID and 2D peaked functions (Gaus- 70 
sians) 0. Vi/2Tin = p Q v 2 {2D} 2 = A Q {2D} defines 
20 quadrupole measure i>2{21?} as a variant of conventional 
v 2 measures which would coincide approximately with 
«2{2} in the absence of jets but which continues to mea- 
sure a nonjet quadrupole (not associated with the SS 2D 75 
peak) to the statistical limits of data in the presence of 
25 jets [13]. The 2D correlation model is not motivated by 
physical interpretations of correlation structure elements. 
The detailed fit model is defined explicitly in Sec. IIIIBI 



refer instead to a nonjet quadrupole (what might be at- 
tributed to elliptic flow) and a jet-related quadrupole (v 2 
contribution mainly from jets and mainly from a promi- 
nent SS 2D peak in 2D angular correlations). 

It has been demonstrated that 2D model fits to angular 
correlations distinguish jet-related structure from the NJ 
quadrupole with few-percent accuracy @, EH, E3, [H, [3(| ■ 
As noted, the model functions used in the 2D model fits 
are motivated by empirical observations of data struc- 
ture, not physical interpretations of structure compo- 
nents 0"S EH, EH- The accurate separation of NJ 
quadrupole and jet-related components by means of 2D 
model fits and estimates of jet- related bias ( "nonflow" ) 
in published v 2 measurements obtained with conventional 
"flow" methods are discussed further in App. [A] 



E. Centrality measures 

Several centrality measures can be related to observed 
charge multiplicity n c h within some angular acceptance 
(A?7, A (ft) based on the Glauber model of A-A colli- 
sions [3l|. The common element is the fractional cross 
section aja^ = b 2 jb\ inferred experimentally from the 
measured minimum-bias event distribution on n c h- 

Glauber Monte Carlo parameters N part (number of 
participant projectile nucleons N) and Nf,i n (N-N binary 
collisions) (with 200 GeV N-N cross section o~nn = 42 
mb assumed) are related to n c h within the STAR TPC 
acceptance. Centrality measure v = 2 Nbi n /N part in turn 
estimates the mean number of N-N binary collisions per 
participant pair (mean participant path length). We re- 
tain the same 200 GeV Glauber parameters for all ener- 
gies as purely geometrical measures of A-A centrality. 



D. "Flow" and "nonflow" 

Conjectured elliptic flow, a possible hydrodynamic re- 

30 sponse to initial pressure/density gradients and overlap 
geometry in non-central A-A collisions, is assumed to be 
measured by the quadrupole (m = 2) term in a Fourier- 
series decomposition of the entire final-state azimuth dis- 85 
tribution [6(. The m — 2 Fourier coefficient is commonly 

35 represented by symbol v 2 . Fourier analysis is applied to 
all azimuth structure ( "anisotropy" ) assuming that el- 
liptic flow dominates that structure. However, possible 
"nonflow" (non-hydrodynamic) contributions to v 2 are 90 
admitted. A variety of schemes has been introduced to 

40 detect and reduce "nonflow" bias. However, considerable 
uncertainty remains for conventional v 2 methods [15| . 
Distinctions between "flow" and "nonflow" have been ex- 
tensively discussed (e.g., Refs. [HI, EE HI])- 95 
The assumptions that support such terminology can 

45 be questioned. Minimum-bias jets contribute strongly to 
"azimuthal anisotropy," are predicted for high-energy nu- 
clear collisions and must form a substantial contribution 
to "nonflow." And the nonjet contribution to v 2 may notioo 
in fact be a flow phenomenon. In the present context we 



F. Two-component model 

The two-component model (TCM) of hadron produc- 
tion near mid-rapidity in high-energy nuclear collisions 
is based on a hypothesis of two dominant production 
mechanisms: (a) projectile nucleon dissociation (soft) 
and (b) large-angle-scattered parton fragmentation to 
jets (hard) [321 ]. The TCM is observed to provide a con- 
sistent quantitative description of pt spectra and corre- 
lations 0, EB~EH HI! • I n more-peripheral A-A collisions 
exhibiting Glauber linear superposition (GLS) of N-N col- 
lisions the soft component should scale cx N part , and the 
hard component (dijet production) should scale oc Nbi n . 

As noted, pair-density difference Ap represents a co- 
variance histogram. If the covariance does not change 
with centrality the ratio Ap/ y/p re f then scales oc l/n c h. 
However, Ap may include contributions from several 
mechanisms with their own scaling behaviors. The soft 
component of Ap/ ^Pref should scale as N part /n c h and 
the hard component as N{,i n /n c h — v x N part /2n c h- That 
is just the scaling observed within the GLS centrality re- 
gion of Au-Au collisions for soft and hard components of 
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2D angular correlations measured by Ap/^p re f Q. 

We can test the TCM more precisely by rescaling the 
fit-model amplitudes used to describe Ap/^p re f with the 
appropriate factors, either n c h/N part (soft) or n c h/Nu n 
(hard). If the rescaled data are invariant on centrality 
over some interval we conclude that the GLS condition 
does persist there and that the assumed soft- or hard- 
component scaling designation is correct. Rescaling re- 
sults are shown in Sec. W\ 



III. 2D ANGULAR AUTOCORRELATIONS 

We describe data volumes, example 2D data his- 
tograms, principal features of the angular correlations 
and fitting procedures used to derive jet-related and non- 
jet quadrupole energy and centrality systematics. 



A. Data histograms 

The analyses reported in Refs. @, were based on 
6.7M and 1.2M Au-Au collisions at y/sj^ = 62.4 (year 
2004) and 200 GeV (year 2001) respectively, observed 
with the STAR TPC. The momentum acceptance was 
defined by transverse momentum p t > 0.15 GeV/c, pseu- 
dorapidity \tj\ < 1 and 2tt azimuth. Au-Au collision cen- 
trality was defined as in Ref. (3lj . Minimum-bias event 
samples were divided into 11 centrality bins: nine ~ 10% 
bins from 100% to 10%, the last 10% divided into two 5% 
bins. The corrected centrality of each bin as modified by ! 
tracking and event vertex inefficiencies was determined 
with a running-integral procedure. Centralitics from N-N 
collisions to central Au-Au were thereby determined to 
about 2% uncertainty. 

Figure [T] (left panels) shows 200 GeV 2D angular corre- 
lations for (corrected) 83-93% (« N-N collisions) and 0- 
5% centrality bins. Angular correlations for 62 GeV have ' 
similar features but with quantitative differences. Within 
the STAR TPC acceptance the minimum-bias correla- 
tion data from Au-Au collisions include three principal 
components: (a) a same-side (SS) 2D peak at the ori- 
gin on (?7a, 4>a) well approximated by a 2D Gaussian for 
all minimum-bias data, (b) an away-side (AS) ID peak 
on azimuth or "ridge" well approximated by AS azimuth 
dipole [1 — cos(0a)]/2 for all minimum-bias data and uni- 
form to a few percent on 77A (having negligible curvature), 
and (c) an azimuth quadrupole cos(20a) also uniform on 
?7a to a few percent over the full angular acceptance of < 
the STAR TPC. Other components consist of a sharp 2D 
exponential peak at (0,0) and a narrow ID peak on r)/\. 
That phenomenological description does not rely on any 
physical interpretation of the components. 

Based on comparisons of observed data systematics ; 
with theory the components (a) and (b) together have 
been in terp reted to represent minimum-bias jets or mini- 
jets fft, Il7{. Component (c), identified as the nonjet az- 
imuth quadrupole, has been conventionally attributed to 




FIG. 1: (Color online) Left: 2D angular autocorrelations for 
200 GeV Au-Au collisions and (a) 83-93% centrality (~N-N 
collisions) and (c) 0-5% centrality. Histograms from 62 GeV 
collisions have the same general features with quantitative 
differences. Right: Two-dimensional model fits to histograms 
in the left panels without the BEC-electron component. 



elliptic flow [FdJ. However, alternative mechanisms have 
been proposed (13 , [3~i| . The narrow exponential repre- 
sents Bose-Einstein correlations and electron pairs from 
photoconversions, and the narrow ID peak on r)& is at- 
tributed to projectile nucleon dissociation. Reinterpre- 
tation of the NJ quadrupole based on comparison with 
jet-related systematics is a main subject of this study. 



Two-dimensional fit model 



In this study we emphasize correlation components (a) , 
(b) and (c). The corresponding 2D model function is p], 



\/Pref 



= An 



Md expi - ~ 



VA_ 

a, 



A D {1 + COS(0A 



®<t> A / \ u 7) A . 

7r)}/2+A Q 2cos(2 A )-(2) 



A ID Gaussian on 77 a (soft component, negligible in 
more-central Au-Au collisions) and 2D exponential (very 
narrow in more-central Au-Au collisions) are omitted 
from Eq. ([2]) for simplicity but were included in the anal- 
yses of Refs. 0, [H I . Equation @ is a more-detailed 
version of Eq. (fj]) . 

Nonjet quadrupole measure Aq as defined by Eq. (T2]) 
is statistically compatible with jet-related measures A2D 
and A]j (all are per-particle measures), permitting 
quantitative comparisons between jet-related and non- 
jet quadrupole systematics. The quadrupole amplitude 
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is related to conventional measure V2{2D} by Aq{2D} = 
Po{b) W 2{2D}, where po(b) — dn cb /2ndr] is the single- 
particle 2D angular density and symbol {2D} denotes 
parameters inferred from 2D model fits to angular corre- 
5 lations as described in Refs. @, H3, ED, HI]. 

IV. TWO-DIMENSIONAL MODEL FITS 

Figure [T] (right panels) shows typical 2D model fits 
compared to corresponding data histograms in the left 
panels. For each data histogram fits are initiated from 

io many different combinations of initial starting parame- 
ters (typically 100-1000) to insure achievement of global 
X minima. The fit residuals are typically consistent 
with bin-wise statistical uncertainties. The general evo- 
lution with centrality is monotonic increase of the SS 2D 

15 peak and AS dipole amplitudes (jet- related structures), 
substantial increase of the SS peak width, rapid de- 
crease to zero of the ID Gaussian on rj^ [7n9| an d non- 
monotonic variation of the nonjet quadrupole 23] . 

A. Jet-related structures 

20 Figure [2] shows fit results for jet-related structures from 
Ref. [3] where they are extensively discussed. The dashed 
curves in the upper panels indicate a Glauber linear su- 
perposition (GLS) trend expected for transparent A-A 
collisions. The jet-related amplitudes follow that trend 

25 from N-N collisions to a sharp transition at v k, 3 cor- 
responding to a /(To ~ 50%. Above that point the am- 
plitudes increase relative to the GLS trend in a manner 
consistent with a modification of parton fragmentation 
that conserves parton energy within resolved jets (l7ll35j. 

30 Figure [2] (lower panels) show the r] and widths of 
the SS 2D peak. Strong elongation on 77 of the SS peak 
in more-central Au-Au collisions was first reported in 
Ref. [8|. The physical mechanism for elongation is cur- 
rently intensely debated [36|, H3] ■ It is notable that the 

35 SS peak azimuth width actually decreases with increasing 
centrality. Mechanisms for jet modification and/or par- 
ton energy loss that rely on multiple scattering and/or 
gluon brcmsstrahlung must confront that decrease. 

B. Nonjet quadrupole 

40 Figure |3] summarizes 2D fit results for Aq{2D}(b) 
(left panel) and corresponding values of V2{2D}(b) (right 55 
panel) for comparison with published V2 measurements. 
The left panel shows fit results for 200 GeV (solid dots) 
and 62 GeV (open circles) data: strong increase to mid- 
45 central collisions followed by reduction to zero for central 
collisions. In the right panel the Vi trend suggests sub- eo 
stantial "elliptic flow" for the most-peripheral centrality 
bin approximating N-N collisions. The solid and dashed 
curves are defined in the next subsection. The dashed 




FIG. 2: Centrality dependence of fit parameters from Eq. © 
for (a) same-side (SS) 2D peak amplitude, (b) away-side (AS) 
ID peak amplitude, (c) SS peak 77 width, (d) SS peak cj> width. 



curves correspond to NA49 «2{EP} measurements at 17 
GeV (see Sec. IIVD|) that provide a reference for inferred 
energy-dependence systematics. 




FIG. 3: Left: Nonjet quadrupole amplitude Ag{2D}(b) plot- 
ted on relative impact parameter 6/60 = \/a/ao with bo ~ 
14.7 fm. Right: Corresponding values for V2{2D}(b). The 
solid and dashed curves are defined by Eq. Q 

The plotting format in the left panel reveals an inter- 
esting centrality trend common to all energies between 
13.5 GeV and 200 GeV: Aq(o) data approximate a Gaus- 
sian trend on relative impact parameter b/bo (&o ~ 14.7 
fm for Au-Au collisions). The dotted curves just visible 
behind the solid curves are modified Gaussians symmet- 
ric about the midpoint 56]. 

The plotting format in the right panel is the more con- 
ventional V2Q)) vs N part (b). Both the correlation mea- 
sure and centrality measure are questionable. Correla- 
tion measure V2 is the square root of per-pair measure v\ 
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actually extracted from pair ratio Ap/p re f. If the physics 
of central Au-Au collisions were not different from p-p 
collisions (Glauber linear superposition) V2 must scale as 
1/ 1 \J n ch- Thus, a general decrease of v% with centrality 
conveys no new information. It is variations relative to 
the 1 / i/rich trend that may be interesting. 

Centrality measure N part biases the visual presenta- 
tion to favor the more-central 50% of the Au-Au frac- 
tional cross section cr/cro- N part = 50 approximates the 
50% point on cr/cr . As shown below, interesting changes 
in jet systematics occur near N part — 50 Q that are ef- 
fectively concealed by the N part centrality measure. The 
lower half of the fractional cross section provides an es- 
sential reference and should be made visually accessible. 
Such access is necessary for proper interpretation of non- 
jet quadrupole data in relation to jet data. Alternative 
centrality measures include b/bo (left panel), fractional 
cross section cr/uo and mean participant path length v. 



is 




0.05- 



^s NN (GeV) 



R(Vs NN ) N hin (b) 



FIG. 4: Left: The energy dependence of azimuth quadrupole 
amplitude Aq — p$v\ evaluated at the maximum of that pa- 
rameter on centrality. The trend 0.13J?(-^/sjvjv) (solid line) de- 
fined above is included in Eq. J3). Right: Nonjet quadrupole 
data (points) for 62 and 200 GeV Au-Au compared to the 
trend defined by Eq. ((3]) (straight line). 



C. Combining energy and centrality trends 



tO yfsNN ' 

in Fig. El 



A3 GeV. It defines the solid and dashed curves 



20 The data in Fig. [3] (left panel) reveal two interesting 
features: (a) Data for all energies above 13 GeV are de- 
scribed by the same centrality variation (solid and dashed 
curves), and (b) the quadrupole amplitude scales with 
energy approximately as log(^/sjviv)- A similar energy 

25 scaling was observed for (p t ) fluctuations/correlations at- 
tributed to minijets [l(J. An algebraic model describ- 60 
ing quadrupole energy and centrality trends (solid and 
dashed curves) is now derived. 

Figure [4] (left panel) shows the trend on collision en- 

30 ergy of Aq{b) maximum values at b/bg ~ 0.5. Published 
i>2 measurements have been converted to Aq values based 65 
on corresponding multiplicity densities and the relation 
Aq = Pqv\. We observe two energy regimes: Below 13.5 
GeV (Bevalac-AGS) the Aq values are small and evolve 

35 with energy from negative to positive in response to the 
kinematic influence of spectator nucleons. The dashed 
line is defined by energy scaling factor R'(-Jsnn) = 
ln(vsjvjv/3.2 GeV) with coefficient 0.008. Above 13.5 
GeV (SPS-RHIC) the rate of increase becomes dramat- 

40 ically larger. The solid line is defined by scaling factor 
Riy/sWW) = \n{y/IWW/!3.5 GeV}/ln(200/13.5) with in- 
tercept 13.5 ± 0.5 GeV and coefficient 0.13 ± 0.01. 

Figure @] (right panel) shows Aq data plotted in the 
form (l/e 2 opt )A Q {2D}(b) vs R(^/s^) N bin (b) , where 

45 Nbin is the number of binary N-N collisions and e op t is 
the eccentricity derived from an optical Glauber simula- 
tion [TH, |38| . We observe empirically that for all Au-Au 
collisions above 13 GeV the Aq{2D} data are described 
accurately by the relation (solid line) 

A Q {2D}{b,^IW^) = C Q R{^I^f) N bm (b) e 2 opt (b).{3) 

so Coefficient C is defined by 1000C = 4.5 ± 0.2. Equa- 
tion © accurately describes the measured p t -integrated 
nonjet azimuth quadrupole in Au-Au collisions for all 
centralities down to N-N collisions and all energies down 



D. Comparisons with other methods 

NGNM V2 measurements are conventionally inter- 
preted to represent some combination of "elliptic flow" 
and "nonflow" [6| with several proposed sources for the 
latter such as resonances. Bose-Einstein correlations and 
jets [Hj]. Separating sinusoids attributed to "flow" from 
other correlation structure is a long-standing problem not 
resolved by non-graphical methods. Various strategies 
have been proposed to reduce "nonflow," including cuts 
on r\ to exclude an interval on t]a containing the origin. 
The most common methods and their biases have been 
compared to 2D model fits on (r? A , 4>a) [II EI [II SI, [11 . 




FIG. 5: Left: Quadrupole results from 2D model fits pre- 
sented in Fig. (solid and dashed curves) compared to pub- 
lished V2 measurements obtained with other methods (points) . 
The dotted curve includes a jet-related contribution obtained 
from the measured SS 2D peak systematics in Fig. [2] Right: 
2D model-fit trend from Eq. Q (solid line) compared to pub- 
lished V2 measurements from other methods (points). The 
dotted and dashed curves include a calculated contribution 
from the measured SS 2D peak. 

Figure [S] (left panel) shows multiparticle cumulant 
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measurements ^4q{2} (open squares) and Aq{4} (solid 
upright triangles) from 200 GeV Au-Au collisions [l5[ 
and event-plane measurements Aq{EP} (solid inverted 55 
triangles) from 17 GeV Pb-Pb collisions [40| compared 
with Eq. ([3]) (solid and dashed curves). In Ref. [14| it 
was demonstrated that W2{2} V2{EP} to within 5% 
(consistent with experiment [15]). "Nonflow" contribu- 
tions to those methods are discussed in App. [A] Pub- eo 
lished uncertainties for the Aq{2} measurements multi- 
plied by factor 10 to make them visible are shown as bars 
within the open squares. It is claimed that «2{4} elim- 
inates "nonnow" arising from small clusters of particles 
(e.g., jets) HU, and those measurements are indeed as 
closer to the Aq{2D} trend (solid curve), but there are 
still very significant differences. The NA49 ^2{EP} mea- 
surements (inverted triangles) (40| provide a reference for 
the energy-dependence systematics. 

Figure [5] (right panel) shows the Aq{2} and Aq{EP} 
measurements compared with Eq. (solid line). Devi- 
ations of the NGNM measurements from the linear trend 
are consistent with expected bias contributions from jet 
structure (SS 2D peak) [l3| as described in App. [AJ The 
dotted curve passing through the 200 GeV measurements 
is a combination of Eq. @ with jet-related contribu- 
tion Aq{SS} derived from jet- related correlation prop- 
erties presented in Sec. IIV Al The same procedure gen- 
erated the dotted curve in the left panel. The dashed 
curve approximating the 17 GeV measurements includes 
the same jet-related Aq{SS} contribution scaled down 
with energy according to R(^snn) from Eq. (|3]). The 
dash-dotted curve indicates what the Aq{2D} data trend 
(solid line) would be if e opt were replaced by Monte Carlo 
ejVfC- The Monte Carlo eccentricity trend appears to 
compensate partially for the jet contribution to Aq{2}. 



show the position of the "sharp transition" (ST) near 
50% (y ps 3 or N part ~ 50). The per-particle peak ampli- 
tudes Ax are rescaled by factor n c h/Nbi n because those 
hard-component amplitudes are expected to scale with 
Nbin (as described in Sec. Ill Fp . The covariance in the 
numerator of Ap/^p re f is then compared directly with 
the number of initial-state N-N binary collisions rather 
than the number of final-state hadrons. We observe that 
below the ST the "jet-related" amplitudes are system- 
atically consistent with a constant value as expected for 
dijet production in a transparent system, thereby but- 
tressing the jet interpretation. Above the ST the ampli- 
tudes increase substantially relative to the GLS trend. 




GLS 
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FIG. 6: Left: Same-side 2D peak amplitude vs fractional cross 
section. The amplitude has been rescaled to compare the 
covariance to the number of binary N-N collisions. Right: 
Away-side ID peak amplitude rescaled in the same way. The 
constant trends are consistent with pQCD dijet production in 
transparent A-A collisions. The hatched bands mark a sharp 
transition in jet-related systematics. 



35 V. QUADRUPOLE VS JET TRENDS 

A unique finding of Ref. Q was the "sharp transition" 
separating Glauber linear superposition and apparent 
A-A transparency within the more-peripheral half of the 
total cross section and strong deviations from the GLS 

40 trend still consistent with a pQCD description within 
the more-central half. The NJ quadrupole data reported 
in Ref. (23j demonstrated remarkably simple energy and 
centrality trends for all Au-Au (or Pb-Pb) centralities 
and all energies above 13 GeV. We make direct com- 

45 parisons between jet-related systematics and the nonjct 
quadrupole. We emphasize the relation of data trends 
to initial-state geometry parameters and possible QCD 
mechanisms. 



A. Minimum-bias jet systematics 

50 Figure |5] shows jet-related SS 2D peak Aid (Gaussian) 85 
and AS ID peak Arj (dipole) amplitudes vs centrality 
measured by fractional cross section. The hatched bands 



In the left panel the SS peak amplitudes for 62 GeV 
rescaled by factor 1.6 [sa 1 / i?(62 GeV) = 1.75] are then 
equivalent to the 200 GeV amplitudes. In the right panel 
the unrescaled AS peak amplitudes for the two energies 
agree closely within the GLS interval. For in-vacuum di- 
jets we expect the SS peak amplitude (sum of individual 
jets projected onto ID y z ) to exhibit a log(^/ s^n) trend 
due to the increase of the kinematically-allowed longitu- 
dinal rapidity y z interval, whereas the AS peak amplitude 
[representing the dijet density on 2D (y z i,y Z 2)] should 
increase more slowly or not at all with energy [3^, |4l[ . 
We observe such trends within the transparency inter- 
val. Thus, comparison of SS and AS centrality and en- 
ergy trends strongly buttresses a dijet interpretation for 
those correlation structures, but also reveals a significant 
change in some jet-related correlation properties above 
the ST. The substantial increase of jet-related amplitudes 
(and SS r\ width) above the ST corresponds quantita- 
tively to possible changes in parton fragmentation that 
still conserve the full parton energy within the resolved 
jet structure 1431 ■ That description is supported by spec- 
trum analysis [l7l [33| and correlation analysis [35| . 
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B. Nonjet quadrupole systematics 

35 

Figure [7] (left panel) shows quadrupole amplitude 
Aq{2D} with the 62 GeV data rescaled by factor 
l/R{&2 GeV) = ln(200/13.5)/ln(62/13.5) = 1.75. The 
close overall agreement is consistent with Eq. ([3]). The 
point-to-point agreement demonstrates the accuracy of 40 
the analysis method, with deviations at the few-percent 
level for two distinct data volumes. The hatched band 
represents the sharp transition in jet properties. It is re- 
markable that in the more-peripheral centrality interval, 
where 3 GeV partons manifest as in- vacuum jets (hadron 
(pt) ~ 1 GeV/c) with no modification and we describe 
A-A collisions as transparent, the nonjet quadrupole con- 
ventionally interpreted to represent "elliptic flow" of a 45 
dense, strongly-interacting QGP increases to 60% of the 
maximum value as measured by Aq{2D} or the maxi- 
mum value as measured by V2{2D} in Fig. [3] (right). 




1 - o/o n 1 - 0/0, 



FIG. 7: Left: Azimuth quadrupole amplitudes vs fractional 
cross section. The 62 GeV data are scaled up according to fac- 
tor 1/R(^/snn) = 1-75 and coincide at the percent level with 60 
the 200 GeV data. The curve is Eq. ©. The hatched band 
indicates that the nonjet quadrupole achieves 60% of its maxi- 
mum amplitude within transparent A-A collisions. Right: Aq 
data rescaled to participant nucleon pairs and further rescaled 
according to Eq. ([3j. Both the rescaled and unrescaled data 65 
are systematically consistent with Eq. ([3|. The dash-dotted 
curve shows the data trend resulting from rescaling with emc 
instead of e opt . 

Fi gure [7] (right panel) shows a variant of Fig. ^ (right 70 
panel) in which the Aq data are first rescaled by factor 
2n c h/N par t to compare the quadrupole covariance in the 
numerator to initial-state participant pairs rather than 
final-state hadrons and then by factor 1/Ni,i n e^ pt corre- 
sponding to Eq. dU). The data are systematically consis- 75 
tent with a constant value (dashed line) . The solid curve 
represents the same rescaling applied to the Eq. ([3]) ex- 
pression for Aq and is also consistent with the data. 

Figure [3] (right panel) shows v 2 {2D} = 0.022 for the 
most peripheral 200 GeV centrality bin (approximately so 
N-N collisions). In Fig. [7] that N-N value is consistent 
with a simple scaling trend describing all Au-Au central- 
ities. It is also notable that V2{2D} in the N-N limit of 
Au-Au collisions is consistent with a pQCD color-dipole 
prediction V2 ~ 0.02 for pions from p-p collisions [34]. 35 



Collision eccentricity can also be modeled by a 
participant-nucleon or Monte Carlo Glauber simula- 
tion [13| . The so-called Monte Carlo eccentricity cmc 
rises well above e op t for peripheral and central collisions 
because of point- like sampling of the nuclear volume flij . 
In Fig. [7] (right) if emc replaced e op t sharp downturns in 
the data would appear at the centrality extremes (dash- 
dotted curve) contradicting any possibility of a conjec- 
tured "hydro limit" for V2 in central Au-Au collisions [45| . 

VI. DISCUSSION 

We summarize and compare several aspects of jet- 
related correlations and jet and nonjet quadrupole data 
systematics and consider alternative interpretations of 
the nonjet quadrupole in light of inconsistencies in the 
conventional "elliptic flow" hydrodynamic interpretation. 

A. 2D model fits compared to other V2 methods 

The azimuth quadrupole Aq{2D} or V2{2D} data used 
in this study are obtained from 2D model fits to angular 
correlations, but the majority of published V2 measure- 
ments are obtained from various nongraphical numerical 
methods (e.g., «2{2}, V2{4}, V2{EP}). It could be ar- 
gued that the 2D fit model is somehow arbitrary, that it 
"depends on assumptions" (e.g., choice of model function 
for the SS 2D peak), and that the inferred quadrupole 
amplitudes are then uninterpretable and should be ig- 
nored [46] . But such arguments imply that NGNM aren't 
based on model fits and don't "depend on assumptions." 

Detailed study of NGNM reveals direct algebraic con- 
nections between such methods and 2D angular corre- 
lation histogra ms, and therefore 2D model fits to those 
histograms [13|, [lj, [3(| ■ In fact, most NGN VI are actually 
based on cosine model fits to 2D angular correlations pro- 
jected onto ID azimuth while subjected to various con- 
ditions on accepted particle pairs, including constraints 
on T) difference acceptance and charge combination. 

The NGNM fit model is effectively a single cosine 
which cannot accurately describe the ID projection. The 
fit residuals are not examined to test the fit validity, and 
the procedure abandons critical information contained in 
the unprojected 2D angular correlations. The fitted co- 
sine amplitude (assumed to represent "flow") can then 
include contributions from multiple correlation compo- 
nents, some identified as "nonflow." Assumptions sup- 
porting NGNM methods include identification of any co- 
sine term as representing a flow [22], [29|, |U |H[ and that 
"azimuthal anisotropy" (any nonuniform structure on az- 
imuth) is dominated by, if not exclusively, flows @ . 

In contrast, the same 2D fit model defined in Refs. @, 
is constrained to describe all 2D data from p-p col- 
lisions to central Au-Au collisions. The model is based 
only on observed prominent features of the data, not on 
physical assumptions. In p-p and more-peripheral Au-Au 
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collisions the SS 2D peak is fully resolved and the data 
require a SS 2D Gaussian model. The remaining struc- 
ture (aside from the soft component and BEC) is fully 
described by two terms of a Fourier series according to in- 

5 spection of the fit residuals. The two Fourier amplitudes 
have very different systematic variations on energy and 
centrality, suggesting minimal parameter covariance. In 
more-central Au-Au collisions the SS peak persists as a 60 
narrow structure on azimuth consistent with a Gaussian. 

10 Systematic uncertainties in the azimuth quadrupole 
arising from 2D model choices are negligible in the GLS 
region where the SS 2D peak is fully resolved. In more- 
central Au-Au collisions it can be shown that the non- 65 
jet quadrupole amplitude is insensitive to the SS peak 

15 r\ structure as long as a SS Gaussian on azimuth is in- 
cluded in the model. Accurate separation of the three ma- 
jor correlation components is confirmed by the internal 
consistency of the parameter trends. Model comparisons 70 
are discussed further in App. [A] 



20 B. Jet correlation systematics 

75 

Figure [5] shows data for the SS 2D peak (left) and AS 
ID peak (right) amplitudes scaled by factor (n c h / Nbi n ) 
to determine the ratio of nominally jet-related covari- 
ances to number of initial-state N-N binary collisions 

25 rather than final-state hadrons. The data are plotted 80 
vs fractional cross section to emphasize an important 
point. Within the lower 50% of the total cross section 
the SS and AS data agree precisely with binary-collision 
scaling as expected for dijet production in transpar- 

30 ent Au-Au collisions, supporting the dijet interpretation. 
The most-probable jets emerge from the lowest-energy 
partons that can appear as jets in the final state (approxi- 
mately 3 GeV), as demonstrated in Refs. 0,[l6rllll- Such as 
low-energy partons (serving as "Brownian probes" Q) 

35 should be most susceptible to a dense, stongly-interacting 
medium. The data are consistent with no jet modifica- 
tion or medium formation in the more-peripheral half of 
the total cross section. * 
Just above the 50% point ( "sharp transition" ) the jet- 

40 related amplitudes increase substantially relative to the 
constant GLS trend, the deviation described as "anoma- 
lous centrality variation" in Ref. [7] . But the increase re- 
mains consistent with pQCD calculations incorporating 95 
modification of fragmentation functions in more-central 

« Au-Au collisions that conserves the parton energy within 
resolved jets [13, The fragment yield increase at 

lower p t (e.g., 0.5 GeV/c) is precisely anti correlated 
with so-called "jet suppression" at larger p t (e.g., 10 lm 
GeV/c) [H|. Jet modification in more-central collisions 

so is not suppression of jet number but rather redistribution 
of fragment number along the jet axis from higher p t to 
lower p t [l3|- We conclude that some aspects of parton 
fragmentation to minimum-bias jets inferred from spec- 
trum analysis and 2D model fits to p t -integral angular 

55 correlations remain consistent with a pQCD jet descrip- 



tion from p-p to central Au-Au collisions. 

C. Quadrupole correlation systematics 

In Fig. [7] (left panel) we demonstrate precise consis- 
tency of 62 and 200 GeV Aq{2D} data scaled by a com- 
mon log(y / SArjv) energy dependence shown in Fig. 0] (left 
panel) and similar to that observed for dijets [lfl]. Com- 
pared to the energy dependence below 13.5 GeV the rate 
of increase above 13.5 GeV is very large (slope changes 
by more than a factor 20). The actual increase in collec- 
tivity below 13.5 GeV is smaller than what the data there 
suggest due to the kinematic effect of spectator nucleons 
( "sqeezeout") resulting in negative values at lower en- 
ergies. Whereas most particles participate in collective 
motion at lower energies, analysis of v-Agt) data to infer 
"quadrupole p t spectra" at 200 GeV [30| suggests that 
only a small fraction of final-state hadrons participates 
in the nonjet quadrupole at higher energies and must 
therefore be very strongly correlated to produce observed 
A Q {2D}(pt,b) amplitudes at RHIC energies [Hill]. 

In Ref. [23| it was demonstrated that Aq{2D} data 
vary approximately as Nbin^-opt- I n the right panel Aq 
is rescaled as {n c h/N part )(l/Nbi n el pt ). The rescaled data 
are again consistent with a constant value (dashed line) 
within ±10% (hatched band), and the data for two en- 
ergies are consistent within a few percent, although the 
absolute Aq values vary over nearly three decades. 

D. Jet-quadrupole comparisons 

By comparing Fig. [S] and Fig. [7] (left panel) we observe 
that Aq attains 60% of its maximum value within a cen- 
trality interval (more-peripheral 50% of cr/<7o) that is ef- 
fectively transparent to jet formation from low-energy 
(mainly 3 GeV) partons, an interval where multiple 
(re)scattering of partons or hadrons apparently plays no 
significant role. Within the transparency interval we ob- 
serve that whereas the jet-related covariance scales as 
Nbin (as expected for pQCD dijet production) the nonjet 
quadrupole covariance scales as N part x Nun x ^tpt > faster 
than dijet production in more-peripheral collisions. 

Figure [7] (right panel) shows that the nonjet 
quadrupole continues to follow the same simple algebraic 
trend within ±10% through and above the sharp transi- 
tion, where the minimum-bias jet trends change dramati- 
cally and where substantial modification of jet formation 
appears [7|. The quadrupole seems to be completely in- 
sensitive to whatever mechanism modifies jet structure. 

E. Implications for hydro interpretations 

What are the implications from these observations for 
hydro interpretations of the azimuth quadrupole? V2 
measurements have been conventionally interpreted in 
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a hydro context in terms of ratio i^/e plotted vs low- 
density limit (LDL) parameter dn c h/dr]S (4f| [5(| (S is 35 
the A-A overlap area). For more-peripheral collisions 
it is expected that vijt oc dn c h/dr]S (assumed corre- 
lated with the mean number of particle rescatterings 
during equilibration). If thermal equilibrium is achieved 
the ideal-hydro limit V2 /e — > constant (saturation) is ex- w 
pected. Previous i>2 measurements were believed to con- 
firm that central Au-Au collisions at 200 GeV achieve 
the ideal-hydro limit (thermalization over some substan- 
tial space-time volume) (5(J HU • 

Conventional vi analysis is based on assumptions that 45 
hydro expansion with particle rescattering is the dom- 
inant dynamical process in heavy ion collisions 0, l5lj . 
that the collision can be described in part as a thermo- 
dynamic state [52|, and V2 is sensitive to an equation of 
state [1, Q . V2 is defined accordingly [H, |(| , and nonflow 50 
contributions to i>2 are estimated with physical-modcl- 
dependent procedures [15j, l2fj, |2l|. The present analy- 
sis presents accurate nonjet quadrupole amplitudes from 
physical-model-independent 2D fits to angular correla- 
tions that reveal simple trends on centrality and collision 55 
energy, including factorization of the dependence on col- 
lision parameters b (impact parameter) and y / S7v r /v- 
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FIG. 8: Left: Ratios V2/e for various combinations of 70 
V2{method} and ex with X — opt or MC. The solid and 
dashed curves are Eq. <(3j for 200 and 17 GeV. «2{2} mea- 
surements for 200 GeV are open squares [15f| [published un- 
certainties (bars) are multiplied by factor 20 for visibil ity] . 
V2{EP} measurements for 17 GeV are solid triangles [4(J. 75 
The dotted lines are Eq. ([3| plus a calculated contribution 
from the jet-related SS 2D peak. Right: Deviations [data 
- Eq. Q] relative to systematic uncertainties for 200 GeV 
2D model- fit measurements (solid dots), and for H2{4} (solid 
triangles), ideal-hydro V2/e = 0.22 (open circles) and LDL 
V2/e = 0.01 (dn c h/drj S) (solid squares). V2{method} values 
are converted to Aq — pqv\. 

Figure |S] (left panel) shows ratio «2/e for several vi w 
methods (v2{2} open squares, V2{EP} solid triangles 
and V2{2D} solid dots) and two e calculations (e op t op- 
tical and €mc Monte Carlo). The vertical scale choice 
excludes a V2{2}/e opt point at 0.55 for central Au-Au 
to retain sufficient resolution for the other data. The as 
data points and line types are consistent with Fig. [5] (left 
panel). The solid curve is Eq. (jHJ including e opt . The 
upper dotted curve is the solid curve plus the contribu- 



tion from the SS 2D (jet) peak calculated in App.lA"! The 
dashed curve represents Eq. © for 17 GeV with e opt , and 
the lower dotted curve is that plus the jet contribution 
from 200 GeV scaled down by factor R(17 GeV). The 
bars in the upper open squares represent the published 
«2{2} systematic uncertainties multiplied by factor 20. 
The lower open squares are the same v 2 {2} measurements 
but combined with tuc- The upper hatched band rep- 
resents an ideal-hydro saturation limit predicted for 200 
GeV. The lower hatched band sketches an LDL trend 
on v. Incorporation of cmc in the W2/e ratio partially 
compensates for the jet contribution to «2{2} and seems 
to meet LDL expectations for more-peripheral collisions. 
But the accompanying downturn for more-central colli- 
sions contradicts hydro expectations for saturation. 

Figure [5] (right panel) compares data and theoretical 
expectations in the form of deviations from Eq. ([3]) di- 
vided by data systematic uncertainties ("error" bars). 
The V2{2D} data (solid dots) are systematically con- 
sistent with Eq. ([3]) as expected. The t>2{4} measure- 
ment deviations relative to their published uncertain- 
ties [5l[ reveal systematic deviations of either sign up 
to ten uncertainty bars from Eq. ([3]). The ideal- hydro 
v 2 /e = 0.22 and LDL-scaling u 2 /e ~ 0.01 (dn c h/dr] S) 
also exhibit large deviations (relative to uncertainties in 
v 2 {2D} data). 

The more-peripheral NJ quadrupole data do not ap- 
pear to follow an LDL trend or require subsequent Au-Au 
collision evolution by particle (parton or hadron) rescat- 
tering, and no transition to an ideal-hydro limiting case 
is observed. The Aq{2D} data do not compel a model 
description based on bulk-medium hydrodynamics or an 
equation of state, in fact contradict such descriptions. 

Thus, the measured Aq{2D} data trends, especially 
the entire energy- dependence trend at and above Be- 
valac energies, the complete insensitivity of the nonjet 
quadrupole to jet modifications in more-central Au-Au 
collisions and detailed understanding of jet-related con- 
tributions to V2 {method} measurements (i.e., required 
corrections for jet-related "nonflow" bias) strongly sug- 
gest that the conventional hydro interpretation of the 
nonjet quadrupole as "elliptic flow" is contradicted by 
most i>2 measurements. 



F. Hydrodynamic flows vs alternative mechanisms 

Other analysis results argue against hydrodynamic 
flows in high energy nuclear collisions |T2| . Published 
V2{pt) measurements for identified hadrons reveal a 
quadrupole p t spectrum common to three hadron species 
and consistent with a boosted source (a form of "radial 
flow") but with a cold spectrum shape very different 
from that for the majority of final-state hadrons [49| . 
The inferred boost distribution is also inconsistent with 
Hubble expansion of a flowing bulk medium [30j|. The 
combined spectrum characteristics suggest that the frac- 
tion of hadrons "carrying" the nonjet quadrupole is sub- 
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stantially less than 10%, ruling out an expanding bulk 
medium as the common source for most hadrons [49j. 
The systematics of v 2 {2D}(p t , b) data reveal that the 
source boost does not depend on Au-Au centrality as 
one might expect for a hydro scenario [30(. < 

Differential study of single-particle p t spectra for iden- 
tified hadrons reveals that spectrum structure conven- 
tionally interpreted with a "blast-wave" model in terms 
of radial flow 53] is actually consistent with parton 
fragmentation to jets for all Au-Au centralities [H, Qjl 
[33|]. The nonjet quadrupole has a substantial ampli- ( 
tude within the "transparency" centrality interval of 
Au-Au collisions, with nonzero values down to N-N col- 
lisions @, [H[ . The measured N-N quadrupole systemat- 
ics extrapolated to LHC energies explain the same-side 
"rid ge" in 7 TeV p-p collisions as a quadrupole manifesta- ; 
tion [4l|,|42]. Thus, the nonjet quadrupole amplitude be- 
comes large in collision systems where particle densities 
are still small, again contradicting a hydro scenario. If 
the nonjet quadrupole does not represent "elliptic flow" 
what is the alternative mechanism? Growing evidence -, 
suggest that the NJ quadrupole is a QCD phenomenon 
arising from low- a; glue-glue interactions leading to long- 
wavelength (multipole) radiation 13j. For example, a 



QCD calculation of interfering radiation from decays 
of two BFKL ladders predicts a long-range quadrupole 
structure in p-p and A- A collisions [54j |. , 
Although the centrality trends for jet-related SS peak 
properties and nonjet quadrupole arc very different the 
two amplitudes, when measured with statistically equiv- 
alent quantities Aq and A 2 u, share the same log( v / SNN) 
energy dependence characteristic of QCD scattering pro- , 
cesses. Equation © reveals that the final-state nonjet 
quadrupole amplitude for (some fraction of) produced 
hadrons is simply determined by initial-state parameters 
(-JsnN) b) over a large kinematic domain including N-N 
(p-p) collisions. There is no evidence for quadrupole sen- , 
sitivity to intermediate processes (multiple rescattering, 
formation of a thermodynamic state, whatever mech- 
anism modifies jet-related correlations above the ST) 
within A-A collisions. And a prediction v 2 ~ 0.02 for 
pions from 200 GeV p-p collisions based on a QCD color- 
dipole model [34| is consistent with the observed N-N ( 
limit v 2 {2D} fa 0.02 in Fig. [3] @, H|. 



VII. SUMMARY 

100 

In summary, 2D angular correlation data for Au-Au 
collisions at 62 and 200 GeV are employed to re-examine 
interpretations of the azimuth quadrupole as "elliptic 
flow" conventionally represented by symbol v 2l a hydro- 
dynamic response to the eccentricity of the initial-state^ 
A-A overlap region. Unlike most conventional v 2 meth- 
ods, two-dimensional (2D) model fits to angular correla- 
tion data distinguish accurately between jet-related fea- 
tures and a nonjet quadrupole represented by symbol Aq. 

The nonjet (NJ) quadrupole exhibits simple sys- 



tematic trends on collision centrality and energy. 
The trend Aq oc fl(v^iViv) N bin (b) e 2 opt {b) with 
R(\/snn) oc log (ysAw/ 13. 5 GeV) accurately describes 
NJ quadrupole data over a broad range of energies and 
all Au-Au centralities. All p t -integrated NJ quadrupole 
data from 17 to 200 GeV are fully described by two A-A 
initial-state parameters. 

In contrast, jet-related features exhibit Glauber lin- 
ear superposition (GLS) trends (transparency) over the 
more-peripheral 50% of the Au-Au total cross section, 
consistent with unmodified dijet production proportional 
to N-N binary collisions. In more-central collisions jet- 
rclatcd amplitudes increase relative to the GLS trend but 
in a manner still consistent with pQCD when a simple 
alteration of fragmentation leading to jets is introduced. 
Within the transparency interval minimum-bias (mainly 
3 GeV) jet characteristics indicate negligible parton or 
hadron rescattering that might contribute to hydrody- 
namic phenomena in a dense medium. Within the same 
centrality interval the nonjet quadrupole amplitude Aq 
increases to 60% of its maximum value. 

Comparison of conventional v 2 measurements with 
jet-related and nonjet-quadrupole systematics reveals 
substantial bias in much of the v 2 measurements due 
to jet-related contributions, the amount depending on 
the v 2 "method" invoked. Comparison of the energy 
and centrality trends inferred for the pt-integral nonjet 
quadrupole with hydro expectations for "elliptic flow" re- 
veals major disagreements. The ratio v 2 /e does not follow 
the number of in-medium rescatterings during equilibra- 
tion (LDL scaling) for more-peripheral A-A collisions and 
does not transition to a near-constant ideal-hydro trend 
for more-central collisions and larger collision energies 
(ideal hydro limit). Quadrupole pt spectra inferred from 
identified-hadron v 2 (p t ) data are very different from the 
spectra for most hadrons. A quadrupole-source boost dis- 
tribution inferred from such measurements is inconsisten 
with Hubble expansion of a flowing bulk medium. 

We conclude that: (i) NJ quadrupole and jet-related 
structures arise from two different mechanisms and can 
be accurately distinguished by 2D model fits, (ii) The 
NJ quadrupole appearing in more-peripheral ( "transpar- 
ent") Au-Au collisions cannot arise from a hydro mech- 
anism relying on many rescatterings. (iii) The same NJ 
quadrupole is insensitive to any mechanism or environ- 
ment (dense medium) that might modify jet structure in 
more-central collisions as observed, (iv) Both jets and 
NJ quadrupole exhibit a similar log( v /s]v : /v) energy de- 
pendence, with threshold near 13 GeV, suggesting a com- 
mon QCD framework for both phenomena, (v) The NJ 
quadrupole is not a hydrodynamic flow manifestation. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DOE under grants DE-FG03-97ER41020 
(UW) and DE-FG02-94ER40845 (UTA). 
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Appendix A: Flows and nonflows 55 

"Elliptic flow" (1)2) measurements rely on an assort- 
ment of analysis methods that encounter the common 
problem of distinguishing "flow" from "nonflow." The 
present study demonstrates that jet structure can be dis- 60 
tinguished accurately from a nonjet remainder that might 
represent "flow" if flow were relevant to high-energy nu- 
clear collisions @, GJ, G3, El HJ- The distinction is 
achieved by 2D model fits to angular correlations. Such 
model fits are dismissed by some because they "depend 65 
on assumptions." But (a) any mathematical analysis de- 
pends on assumptions — the question is the validity of 
what is assumed — and (b) NGNM measurements also 
rely on model fits. 

The NGNM V2 fit model is a cosine (or cosine plus 70 
constant) applied to a ID projection onto azimuth of all 
2D angular correlations. Different V2 methods are dis- 
tinguished by the conditions imposed on accepted pairs 
(charge combination, hadron species, rj acceptance) in 
attempts to reduce "nonflow" bias to v m based on physi- 75 
cal assumptions. The fit residuals are not presented (but 
see comments on ZYAM subtraction below). Based on 
results from Ref. Q the residuals for a ID single-cosine 
model must be large. Some fraction of the jet struc- 
ture must be included in NGNM V2 measurements as a 
"nonflow" bias. The amount of jet-related bias depends 
on the "method." The bias can be predicted accurately 80 
from 2D model fits if V2{method} is sufficiently well de- 
fined @, m . 

1. System A vs System B 

85 

We can identify two descriptive systems. System A 
is based on the observation that within all 2D angular 
correlation data three prominent features or components 
labeled (a), (b) and (c) consist of a SS 2D peak, an AS 
dipole and an azimuth quadrupole not associated with 
(a). Those features persist for all A- A collision systems. 
We measure their characteristics with 2D model fits. No 
physical assumptions motivate that description. 

pQCD provides the standard description of high en- 
ergy nuclear collisions, a falsifiable theory that makes 95 
real predictions about what should be observed. Included 
in those predictions is the appearance of (a) and (b) in 
2D angular correlations and their systematic properties. 
The correspondence between predictions and data re- 
quires interpretation of (a) and (b) as pQCD jet-related. 100 
The jet interpretation arises from and relies on an estab- 
lished physical theory. What remains is (c), the nonjet 
quadrupole distinct from (a) and (b) in all cases. Its in- 
terpretation is questioned. No falsifiable theory currently 
predicts all measured properties of the NJ quadrupole. 
Thus, System A can be represented by Data = pQCD 
"jets" + "nonjets" and is the basis for Ref. Q reportingios 
minijet systematics. 

System B is based on the primary assumption that 



"flows" (collective motion shared by many particles) 
must play a major role in high energy nuclear collisions. 
But "flows" are not required by QCD to exist at higher 
energies. The nucleons participating in collective motion 
observed at the Bevalac are no longer relevant. High- 
energy collisions are dominated by the low-x gluons in 
the projectiles - either liberated in place to form "soft" 
hadrons or undergoing "hard" scattering to form dijets. 

In System B "Flows" are associated with one or more 
cylindrical multipoles. Non-graphical numerical methods 
(NGNM) based on physical assumptions are developed to 
extract multipole amplitudes v m interpreted to represent 
"flows." The methods rely on a common assumption that 
"flows" can be modeled by cosines. NGNM methods (the 
fit models) do not recognized the existence of a SS (jet) 
peak in correlation data or the contributions that the 
SS peak must make to NGNM multipole amplitudes v m . 
But the SS peak is the dominant source of "nonflow" in 
published v m measurements. The network of v m meth- 
ods and assumptions is complex and changeable. Sys- 
tem B is represented by Data = "flow" + "nonflow" + 
"other" ( "other" represents structure independent of the 
v m "flow" multipole of interest) and is the basis for anal- 
ysis in Ref. (46j and other dihadron correlation studies 
based on ZYAM subtraction 39]. 

2. Mapping from System A to System B 

If a "flow" analysis method is sufficiently well-defined 
we can establish a quantitative relation between systems 
A and B and demonstrate that the cross terms are large 
for the usual "flow" methods. "Jets" from A is split be- 
tween "flow," "nonflow" and "other" in B. "Nonjets" is 
also split, and some of "nonjets" may appear in "non- 
flow." There is no justification for assuming that "non- 
flow" -(- "other" includes all of "jets," but that common 
assumption is the basis for ZYAM subtraction [39|, H(| . 

All angular correlation data include a SS 2D peak [fea- 
ture (a), part of "jets"]. The SS peak is always narrow 
on <j) but is elongated on r\ in more-central Au-Au col- 
lisions [7j. The SS 2D peak projected onto ID azimuth 
is actually a periodic peak array with a Fourier series 
equivalent. The Fourier amplitudes for given peak width 
are represented by factor F m (a^, A ) for the m th Fourier 
term (cylindrical multipole) [39|. Projection of the 2D 
peak onto ID azimuth depends on its rj width relative 
to the rj acceptance Ar] and is represented by factor 
G(a v& , Arj) [33. The jet-related quadrupole amplitude 
derived from SS 2D peak properties is then defined by [22[ 

2A Q {SS}(b) = F m (a^)G(a m ,A V )A 2D (b). (Al) 

With that expression we can relate nonjet Aq{2D} and 
jet- related Aq{SS} quadrupole amplitudes in System A 
to "flow" and "nonflow" in System B. The expression for 
more-complex 77-exclusion cuts is derived in Ref. [22j . 

Figure |H] (left panel) shows a parametrization of the 
centrality dependence of SS 2D peak amplitude A2d(1>) 
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and three quadrupole amplitudes related by [l3[ 

A Q {2}=A Q {2D}+A Q {SS}. (A2) 

A q {2D} (solid curve) is denned by Eq. © [H| and 
Aq{SS} (dashed curve) by Eq. (|Al|) using SS peak pa- 
rameters (amplitude and widths) from Ref. [7j. ^4q{2} 
(dotted curve) is then a prediction for measured v 2 {2} ~ 
V2{EP} measurements that are derived from cosine fits 
to ID projection onto azimuth of all 2D angular correla- 
tion structure ("flow") QUI]. The dotted curve in the 
left panel appears (transformed) in the right panel and 
Fig. El 




FIG. 9: Left: SS 2D (jet) peak amplitude Aid (parametriza- 
tion of data from Ref. [3]), SS 2D peak quadrupole compo- 
nent Aq{SS} ("nonflow," inferred from Ref. [7| data) and 
nonjet quadrupole amplitude Aq{2D} from Eq. with 60 
A Q {2} = A Q {2D} +A q {SS} Q. Right: Quadrupole am- 
plitudes Aq{X} converted to conventional measures V2{X}. 
The open squares are i>2{2} measurements from Ref. [l5| . The 
curves correspond to those in the left panel. 

65 

Figure [S] (right panel) shows t>2{X} trends obtained 
from the corresponding Aq{X} curves in the left panel 
by Aq — pof|. Also included are v 2 {2} measurements 
(open squares) from Ref. [l5j . An equivalent compari- 
son for 17 GeV measurements appears in Fig. [5l The ™ 
precise agreement between measurements (points) and 
prediction (dotted curve) is evident. From the nonjet 
quadrupole trend and jet-related correlation structure in 
System A we accurately predict V2{EP} « v 2 {2} pub- 
lished "flow" measurements in System B. The prediction 75 
does not include small contributions to «2{2} from BEC 
and electron pairs that are excluded from the SS peak 
A 2 r>{b) data by the 2D model- fit procedure Q- For sta- 
tistically well-defined v 2 methods (e.g., v 2 {2} w v 2 {EP}) 
the "nonflow" (jet) bias contribution to inferred "flow" so 
v 2 can be estimated accurately. 

Attempts have been made to parameterize "nonflow" 
contributions to v 2 with quantity g 2 = N part [v 2 {2} — 
v|{4}] [l5( | . If we approximate «2{4} by v 2 {2D} we 
obtain 52 = {N part /p )A Q {SS} « {5/[l + 0.l(v - 35 
1)]}Aq{55}. The factor in curly brackets is derived from 
a two-component model of particle production [32] . That 
expression agre es q uantitatively with 200 GeV g 2 data in 
Fig. 31 of Ref. [15| except for more-peripheral collisions 
where Aq{SS} derived from 2D model fits does not in- 90 
elude a relatively large contribution to W2{2} from BEC 



and electron pairs. The comparison confirms the large 
"nonflow" contribution to «2{2} from the SS jet peak. 



3. Consequences of System B 

System B is the basis for dihadron correlation analy- 
sis on ID azimuth including "ZYAM subtraction" of a 
combinatoric background. A "flow" background is esti- 
mated based on NGNM v 2 measurements and subtracted 
from Data ("raw" correlations). It is assumed that for 
some combination of v 2 methods the difference Data — 
"flow" = "nonflow" + "other" includes all of "jets." The 
subtraction does reveal the large residuals of the System 
B model fit. After application of "trigger-associated" pt 
cuts it is assumed that the surviving "nonflow" + "other" 
still retains all of "jets." Since "flow" actually includes 
some fraction of "jets" ("nonflow") the surviving "jets" 
structure in the ZYAM-subtracted and pt-cut "nonflow" 
correlations is attenuated and distorted, leading to incor- 
rect inferences about jet systematics [39(. The "nonflow" 
+ "other" component in System B may include some 
fraction of "jets" from System A, but the fraction de- 
pends on arbitrary definitions of "flow" and "nonflow." 

Reference [46j presents a direct comparison between 
System A and System B. The ZYAM subtraction (Sys- 
tem B) shown in Figs. 3 and 7 seems to indicate that with 
increasing angle relative to the EP from in-plane to out- 
of-plane the jet structure is increasingly attenuated and 
distorted, implying that "jet quenching" is directly cor- 
related with the apparent parton path length in a "dense 
medium." What survives ZYAM subtraction is then fur- 
ther separated into "jet-like" and "ridge-like" structure, 
again based on arbitrary assumptions. The paper con- 
cludes "...high p t triggered jets are biased toward surface 
emission, and the jet fragmentation is hardly modified by 
the medium" (i.e., jets in central Au-Au collisions are the 
same in structure and abundance as in p-p collisions). 

But Fig. 4 of that paper presents results from ID model 
fits (System A), including a SS peak model element, that 
reveal undistorted jet structure. In contrast to severe jet 
attenuation and distortion with increasing angle inferred 
from System B the System A model fits reveal a possible 
increase of jet correlation amplitude and no distortions, 
consistent with the minijet analysis in Ref. Q in which 
minijets are modified in more-central collisions, but in a 
manner consistent with pQCD expectations (l7l l33l. |35| . 
The System A results in Fig. 4 of Ref. [46| are dismissed 
there because v 2 systematics inferred from ID model fits 
contradict "flow" v 2 inferred from NGNM analysis. 

Any application of v 2 measurements to subsequent 
analysis requires a choice between (a) a comprehensive 
2D fit model motivated by actual data structures com- 
bined with a few physical-model-independent assump- 
tions and (b) an incomplete cosine model (NGNM) com- 
bined with many physical-model-dependent a priori as- 
sumptions. Choosing (a) may lead to interesting new 
physics derived within an intact QCD context represent- 
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ing a falsifiable theory. Choosing (b) (e.g., ZYAM sub- leading, 
traction) can lead to results that are distorted and mis- 
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